Higher-twist contributions to large p ± hadron production in hadronic collisions 
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The scaling behavior of large-p^ hadron production in hadronic collisions is investigated. A 
significant deviation from the NLO QCD predictions is reported, especially at high values of 
x ± = 2p ± /y r s. In contrast, the prompt photon and jet production data prove in agreement 
with leading-twist expectations. These results are interpreted as coming from a non-negligible 
contribution of higher-twist processes, where the hadron is produced directly in the hard 
subprocess. Predictions for scaling exponents at RHIC are successfully compared to PHENIX 
preliminary measurements. We suggest to trigger on isolated large-p^ hadron production to 
enhance higher-twist processes, and point that the use of isolated hadrons as a signal for new 
physics at colliders can be affected by the presence of direct hadron production processes. 



1 Introduction 

The most important discriminant of the twist of a perturbative QCD subprocess in a hard 
hadronic collision is the scaling of the inclusive invariant cross section^, 

o*" = E%-{AB 4CX) = ^ , (1) 
dV P" 

at fixed x x = 2p x /y/~s and center-of-mass (CM) angle In the original parton model the power 
fall-off is simply n = 4 since the underlying 2 — > 2 subprocess amplitude for point-like partons 
is scale invariant, and there is no dimensionful parameter as in a conformal theory. However, in 
general additional higher-twist (HT) contributions involving a larger number of elementary fields 
contributing to the hard subprocess, n act i V e > 4, are also expected. For example, the detected 
hadron C can be produced directly in the hard subprocess as in an exclusive reaction. Unlike 
quark or gluon fragmentation, the direct processes do not waste same-side energy, thus involving 
minimal values of the momentum fractions x\ and X2 where parton distributions are maximal. 
Neglecting scaling violations in QCD, the scaling exponent n is given by n = 2n ac ti V e — 4. 
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The idea of direct hadron production was considered in the 1970's to explain the large fixed 
x x scaling exponents reported at ISR and fixed target FNAL energies^. However, there has been 
no comprehensive and quantitative analysis of the data up to now which could bring compelling 
evidence for such higher-twist effects. In these proceedings, we summarize the novel aspects 
discussed in our recent analysis^, namely: 

(i) a dedicated analysis of the most recent FNAL, RHIC and Tevatron data on large-p ± 
hadrons, prompt photons and jets; 

(ii) the systematic comparison of the experimental scaling exponents with NLO QCD expec- 
tations; 

(iii) predictions for the top RHIC energy and at the LHC. 



2 Analysis 

The exponent n NLO of mid-rapidity particle production has been computed in QCD at next-to- 
leading order (NLO) accuracy from Ref.^1. The x x -dependence of n NLO at fixed p ± has been 
determined for various hadron species (ir, K, p/p). At p ± = 10 GeV the exponents increase 
slowly from n NLO ~ 5 at small values of x ± (x ± = 1CP 2 ) up to n NLO ~ 6 at x ± = 0.5 with 
almost no dependence on the specific hadron species. Remarkably, the exponent extracted in the 
prompt photon channel is below those of hadrons, by roughly one unit, close to the conformal 
limit, n = 4, at the smallest values of x ± . This observation is understood from the (relative) 
absence of fragmentation processes and one less power in a s , leading to less scaling violation in 
this channel. 

On the experimental side, the exponent n exp has been systematically extracted from mea- 
surements in hadronic collisions, from fixed-target to collider experiments. It is deduced from 
the comparison of x ± -spectra at different CM energies, 

.xp, > = ln(a inv (x ± ,yir)/a-(^,y^)) 
In \^fsil^Js%) 

which is equivalent to (pQ) at fixed x ± . The data sets include 7r° measurements by the E706 at 
fnalSI and by the PHENIX collaboration at RHIC El. At higher energies, the measurements 
of charged hadrons (or charged tracks) in p-p collisions at yfs = 630, 1800 GeV by CDF^ and 
= 500, 900 GeV by UAl^ are included in the analysis. Also considered are prompt photon^ 
and jet^l data obtained by CDF and DO at = 546, 630, 1800 GeV. 



3 Results 



The hadron exponents plotted in Fig. [T] (left) exhibit a clear trend, with a significant rise of n exp 
as a function of x ± . Typical values of n exp are n exp ~ 5-6 at small x ± ~ 10~ 2 while PHENIX 
data point to a mean value n exp ~ 6.7 at an intermediate x ± ~ 10 _1 . At higher values of x ± , 
the measurements by PHENIX and E706 reveal an exponent even larger, n exp ~ 8, confirming 
the results reported long ago at the ISR. The exponents obtained in the photon and jet channels 
are strikingly different, showing almost no dependence on x ± . Importantly enough, the values 
obtained lie only slightly above the conformal limit, ~ 4.6 and nj* t g ~ 4.4, i.e. several units 
smaller than the exponents observed for hadrons. 

In order to compare properly data and theory, the difference between experimental and 
theoretical exponents, A(x ± ) = n exp — n NLO , is plotted in the right panel of Fig. [T]for hadrons 
and photons/jets. Note that the error bars include both experimental as well as theoretical 




errors, added in quadrature. The theoretical uncertainty is estimated from the variation of 
renormalization/factorization scales fromp x /2 to 2p ± . Fig. [1] (right) indicates that the hadronic 
exponents extracted experimentally prove significantly above the leading-twist (LT) predictions. 
The discrepancy is moderate at small x x ~ 10~ 2 , A ~ 0.5, but becomes increasingly larger at 
higher values of x ± : A ~ 1 at x ± = 10 _1 and up to A ~ 2 in the largest x ± region. In contrast, 
the scaling behavior observed for photons and jets are in very good agreement with the NLO 
predictions (A ~ 0). 



4 Discussion 



Part of the discrepancy reported in hadron production data at large x ± ~ 1 could occur because 
of the appearance of large threshold logarithms, ln(l — x ± ), which should be resummed to all 
orders in perturbation theory^l. It would therefore be most interesting to investigate whether 
or not threshold resummation might bring data and theory in agreement. Note however that 
the discrepancy is also observed at small values of x ± ~ 10 -2 , where such effects are usually 
expected to be small. 

A natural explanation for the large exponents observed in the hadron channel is the presence 
of important HT contributions from processes in which the detected hadron is produced directly 
in the hard subprocess, because of the dimension of the hadron distribution amplitude. In 
contrast, particles having no hadronic structure like isolated photons and jets are much less 
sensitive to such HT contributions and should behave closer to LT expectations, as observed. 
Another piece of evidence for HT effects is the larger exponents for protons than for pions 
observed at the ISR. As discussed in^, the difference between the direct proton and pion scaling 
exponent is n p — = 2 (n p = 8, n n = 6) instead of n p — n n ~ at LT. The experimental 
value obtained from the ISR, n p — ~ 1 , thus reflects the mixture of LT and HT contributions 
to the total cross section. It has also been noted^^ that the presence of color-transparent HT 
subprocesses can account for anomalous features of proton production in heavy ion collisions^! 

Finally, we discuss the phenomenological consequences of possible HT contributions to 
hadron production in p—p collisions at RHIC and LHC. In order to obtain qualitative pre- 
dictions, the difference A between the experimental and the NLO exponent has been fitted 
to the hadron data currently analyzed. The typical values of A fit expected at RHIC (taking 
y/s = 200,500 GeV) and at LHC (y/s = 7 TeV, compared to s/s = 1.8 TeV at Tevatron) are 




plotted as a function of x x in Fig. H At RHIC, A fit is slightly below 1 at small x ± < 5.10" 2 but 
decreases towards zero at larger x ± . The predictions turn out to be in very good agreement, 
both in shape and magnitude, with the PHENIX preliminary measurements ^21 performed at 
y/s = 500 GeV. At LHC, smaller deviations with NLO expectations are expected because of the 
large values of p ± probed at high energy: A fit ~ 0.5 below x ± =5xl0 3 and smaller above. In 
order to enhance the HT contribution to hadron production, we suggest to trigger on isolated 
hadrons, i.e. with small hadronic activity in their vicinity. We also point that the use of isolation 
cuts, usually applied for prompt photons, will strongly suppress LT processes. Consequently, 
the scaling exponents of isolated hadrons are expected to be somewhat larger than those in the 
inclusive channel. The use of isolated hadrons as a signal for Higgs production or new physics 
scenarios^^ might be confused by the presence of direct hadron production processes. 
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